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PREFACE 

The  work  reported  herein  is  part  of  a  basic  research  on  trailing 
vortices  in  stratified  and  unstratified  fluids  and  was  performed  during 
the  period  30  December  1983  to  30  March  1984.  The  principal  investigator 
is  Dr.  Turgut  Sarpkaya,  Distinguished  Professor  of  Mechanical  Engineering 
at  the  Naval  Postgraduate  School.  Sponsorship  was  provided  by  the 
Defense  Advanced  Research  Projects  Agency  (DARPA)  under  ARPA  Order  3925. 
The  contract  technical  monitor  was  Dr.  Arthur  J.  Bruckheim,  Program 
Manager,  Tactical  Technology  Office,  DARPA. 
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I.  INTRODUCTION 


Vortices  and  vortex  wakes  have  become  a  major  theme  of  aerodynamics 
research  since  the  advent  of  the  large  aircraft  and  the  understanding  of 
their  evolution  required  an  examination  of  many  fundamental  problems  in 
fluid  mechanics.  Much  of  the  progress  made  during  the  past  two  decades 
was  discussed  at  the  Symposium  on  Aircraft  Wake  Turbulence  and  Its 
Detection  [Ref.  1]  and  at  the  Aircraft  Wake  Vortices  Conference  [Ref.  2]. 
Comprehensive  reviews  of  the  entire  subject  have  been  given  by  Donaldson 
and  Bilanin  [Ref.  3],  Widnall  [Ref.  4]  and  Hal  lock  and  Eberle  [Ref.  5]. 

These  studies,  as  well  as  numerous  others  carried  out  since  1977, 
have  uncovered  a  number  of  complex  problems  which  must  be  resolved  in 
order  to  achieve  a  better  understanding  of  the  important  features  of 
trailing  vortices  in  homogeneous  and  stratified  media.  The  principal 
ones  are  as  follows  [Ref.  6]: 

(i)  Roll-up  process:  The  velocity  and  turbulence  distributions  at 
any  station  behind  the  wing  depend  on  the  wing  section,  wing-tip  shape, 
Reynolds  number,  wing  incidence,  and  the  distance  of  the  station  from 
the  wing  [Ref.  7].  The  distributions  of  the  initial  velocity  and  turbu¬ 
lence,  which  influence  the  roll-up  and  the  decay  process,  cannot  be  changed 
independently.  For  example,  a  change  in  tip  shape  changes  the  core  size, 
as  well  as  the  velocity  and  turbulence  distributions.  High  levels  of 
turbulence  result  in  an  increased  diffusion  of  vorticity,  which  in  turn 
increases  the  core  size. 


(ii)  Probe  sensitivity  of  the  vortices:  Flow  visualization  studies 
suggest  that  trailing  vortices  are  extremely  sensitive  to  disturbances 
created  by  even  very  small  probes  or  bubbles.  This  forces  one  to  use 
non-intrusive  means  of  measurements  such  as  an  LDV.  Even  then,  'vortex 
wandering'  [Ref.  8],  which  makes  the  vortices  appear  larger  than  normal 
in  time-averaged  velocity  measurements  (for  vortices  generated  by  a  wing 
in  a  wind  tunnel),  or  the  unsteady  nature  of  the  flow  (for  vortices 
generated  by  a  wing  in  a  tow  basin)  makes  the  mean  velocity  profiles  in 
the  vortices  difficult  to  determine. 

(iii)  Large-scale  instabilities:  The  vortices  are  seldom  observed  to 
decay  away  owing  to  viscous  and  turbulent  dissipation,  but  are  almost 
always  destroyed  by  either  mutual  induction  instability  (Crow  instability 
[Ref.  9])  and/or  vortex  breakdown.  The  Crow  instability  grows  exponen¬ 
tially,  and  results  either  in  a  linking  of  the  vortex  pair  into  a  series 
of  crude  vortex  rings  or  in  a  highly  disorganized  intermingling  of  the 
vortices.  Vortex  breakdown,  whose  mathematical  details  have  not  yet 
been  adequately  treated,  rearranges  the  vortex  structure  and  increases 
the  core  size,  turbulence,  and  energy  dissipation.  Thus,  it  is  very 
difficult  to  measure  accurately  the  trajectories  of  the  three-dimensional 
vortices  from  their  creation  to  their  ultimate  demise. 

(iv)  Reynolds  number:  Even  the  highest  Reynolds  numbers,  based  on 
wing  chord,  reached  in  wind  tunnels  or  towing  basins,  are  an  order  of 
magnitude  lower  than  what  is  possible  for  an  aircraft.  Thus,  the  scale 
effects  are  not  easy  to  assess. 

(v)  Ambient  conditions  such  as  turbulence  and  stratification  play 

major  roles  in  the  evolution  of  vortices.  The  quantification  of  these 

effects  offers  exceedingly  complex  difficulties.  Finally, 
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(vi)  Ground  or  free-surface  effects:  The  Vortex  pair  may  move 
towards  a  rigid  boundary  at  which  the  no-slip  condition  must  be  satis¬ 
fied  or  towards  a  free  surface  at  which  the  zero-shear  condition  must  be 
satisfied.  In  either  case,  the  vortices  come  under  the  influence  of 
their  images  and  the  mutual  induction  instability  may  occur  not  only 
between  the  two  real  vortices  but  also  between  a  real  vortex  and  its 
image.  In  general,  one  would  expect  mutual  induction  instability  among 
four  vortices  (two  real  and  two  image  vortices). 

The  phenomenon  is  further  complicated  by  several  additional  facts. 
When  the  vortices  are  propelled  towards  a  rigid  surface,  vorticity  of 
opposite  sign  is  generated  on  the  no-slip  boundary  and  swept  towards  the 
vortex  pair.  The  total  vorticity  diminishes  very  quickly  as  vorticity 
from  the  two  regions  diffuses,  the  wall  region  serving  as  a  strong  sink 
for  the  vorticity  associated  with  the  original  vortex  [Ref.  10].  The 
development  of  a  boundary  layer  along  the  rigid  wall  may  give  rise  to 
flow  separation  for  sufficiently  high  Reynolds  numbers.  With  or  without 
such  a  separation,  however,  the  center  of  the  vortex  pair  eventually 
moves  away,  or  'rebounds',  from  the  wall  [10-12]. 

For  the  case  of  a  zero-stress  boundary,  the  free  surface  still  acts 
as  a  vorticity  sink,  but  this  is  relatively  weak  due  to  the  absence  of 
intense  oppositely-signed  vorticity.  Thus,  in  the  absence  of  other  im¬ 
peding  phenomena,  one  expects  a  mild  interaction  between  the  vortices 
and  the  free  surface  and  a  small  rebound  of  the  vortex  center  from  the 
free  surface.  However,  the  ability  of  the  free  surface  to  deform  under 
the  influence  of  strain  fields  leads  to  a  strong  interaction  between 
the  vortices  and  the  free  surface. 
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When  the  vortex  pair  approaches  the  free  surface  (with  or  without 
mutual  induction  instability  and/or  vortex  breakdown),  the  vortices 
(or  the  crude  vortex  rings)  and  their  images  give  rise  to  distinct 
surface  disturbances  (waves  of  very  small  length).  These  waves  or  the 
local  roughening  of  the  free  surface  are  visible  because  the  particle 
motions  they  create  as  they  propagate  produce  a  straining  field  at  the 
surface  which  modifies  the  physical  characteristics  of  the  sea  surface 
(reflection,  refraction  and  the  scattering  of  sound,  heat  and  light). 

In  fact,  in  many  photographs  of  the  sea  surface  taken  from  aircraft  or 
satellites  using  passive  or  active  sensors,  the  differential  roughening 
of  the  free  surface  can  be  seen  clearly  [Ref.  13].  It  is  important  to 
note  that  the  roughening  of  the  sea  surface  is  brought  about  not  only 
by  the  proximity  of  vortices  or  vortex  rings  but  also  by  internal  waves 
[Ref.  13].  However,  the  character  of  the  surface  signature  strongly 
depends  on  the  manner  it  is  generated  and  there  is  no  difficulty  in 
distinguishing  the  surface  signature  created  by  the  vortices  from  that 
created  by  the  internal  waves. 

The  surface  signatures  produced  by  a  pair  of  trailing  vortices 
generally  occur  in  the  later  stages  of  the  vortex  lifetimes,  and  as  such 
are  subject  to  the  complex  mechanisms  previously  cited.  Thus,  it  is 
essential  to  achieve  some  understanding  of  the  evolution  of  vortices 
in  an  infinite  homogeneous  or  stratified  medium  prior  to  embarking  on 
an  investigation  of  the  free-surface  effects  and  surface  signatures. 

A  number  of  experimental  and  analytical  studies  has  been  carried  out 
at  the  Naval  Postgraduate  School  by  Sarpkaya  and  his  students  [Refs.  6, 
14-16]  in'order  to  Investigate  the  effects  of  the  governing  parameters. 
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with  the  exception  of  those  due  to  surface-proximity,  on  the  rise  and 
demise  of  the  trailing  vortices  in  homogeneous  and  density  stratified 
media.  These  studies  have  clearly  identified  the  various  demise  mech¬ 
anisms  in  both  media  and  established  basic  relationships  between  the 
rise  of  vortices  and  the  governing  parameters  in  an  effectively  infinite 
medium.  It  is  against  this  background  that  the  present  investigation 
of  the  free-surface  effects  and  the  evolution  of  surface  signatures 
is  undertaken. 


A.  EQUIPMENT 

The  equipment  used  to  generate  the  trailing  vortices  has  been  exten¬ 
sively  used  at  this  facility  over  the  past  three  years  [Refs.  6,  14-16]. 
Only  the  salient  features,  most  recent  modifications  and  the  adaptation 
for  this  work  are  briefly  described  in  the  following. 

The  system  consists  essentially  of  a  towing  basin.  The  auxiliary 
components  of  the  basin  are  the  plumbing  for  water  and  for  water  and 
brine  mixture  (to  establish  a  density  stratification),  turbulence  manage¬ 
ment  system,  top  and  bottom  carriages,  velocity  measuring  system,  light¬ 
ing  system,  and  the  models  [Refs.  14-16]. 

A  mirror,  inclined  45  degrees  from  the  horizontal,  is  attached  to  the 
top  of  the  tank  at  the  test  section  to  allow  viewing  of  the  vortex  pair 
and  the  surface  disturbances  from  the  top  of  the  basin.  As  the  shadows 
of  the  surface  disturbances  rather  than  the  disturbances  themselves  were 
photographed  (taking  advantage  of  the  differential  roughening  of  the 
water  surface  and  the  change  of  its  reflective  properties),  a  single 
light  was  positioned  along  the  axis  of  the  basin  at  a  distance  far 
enough  to  reduce  glare  and  parallax  and  near  enough  to  produce  sharp 
shadows  of  the  surface  disturbances.  Sheets  of  plexiglas  (painted  flat 
black)  were  placed  on  the  side  walls  of  the  test  section  to  eliminate 
undesirable  reflections  (see  Fig.  1). 

The  drains  are  provided  at  the  bottom  at  each  end  of  the  basin.  Two 
parallel  rails  are  mounted  along  the  bottom  of  the  tank.  A  carriage 
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rides  smoothly  on  these  rails  and  provides  the  test  body  with  a  constant 
velocity  through  the  use  of  an  endless  cable  and  a  DC  motor.  The  velocity 
of  the  model  is  measured  and  continuously  monitored  through  the  use  of  a 
magnetic  linear  displacement  transducer.  It  yields  a  signal  proportional 
to  the  velocity  of  the  model  with  an  error  less  than  1  percent. 

The  two  rails,  the  carriage  and  the  filling  pipes  are  located  on  or 
near  the  bottom  of  the  basin  and  under  a  turbulence  management  system. 

The  system  consists  of  one  inch  thick  polyurethane  foam  sandwiched 
between  two  perforated  aluminum  plates. 

B.  MODELS 

Three  models  were  used  in  the  experiments:  A  rectangular  foil  (NACA 
0012)  and  two  Delta-wings,  one  with  sharp  edges,  and  one  with  rounded 
edges.  These  are  the  same  models  which  were  used  previously  at  this 
facility  [Refs.  6,  14-16].  The  interior  of  each  model  was  hollowed  and 
used  as  a  dye  reservoir  to  seed  the  vortex  cores.  The  models  are  mounted 
on  their  bases  by  means  of  a  thin  streamlined  aluminum  bar  with  a  cross 
section  of  a  NACA  0006  foil  and  set  at  the  desired  angle  of  attack.  As 
noted  earlier,  all  models  are  pulled  by  means  of  a  DC  motor,  pulley,  and 
cable  system  at  the  desired  speed  {ranging  from  1.0  to  5.0  feet/second). 

C.  TEST  PROCEDURES 

The  model  is  placed  in  the  basin  and  the  basin  is  filled  gradually 
with  fresh  water  to  the  desired  level.  In  the  cases  where  dye  was  used 
to  examine  the  vortex  pair,  the  models  were  filled  with  dye  prior  to 
filling  the  tank.  After  removal  of  trapped  air  and  after  a  sufficient 
period  of  waiting  for  the  escape  of  dissolved  air  in  the  water  and  the 
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elimination  of  any  internal  currents  in  the  basin,  the  model  was  set 
in  motion  at  the  desired  speed. 

The  shadowgraph  technique  was  used  to  determine  the  characteristics 
of  the  surface  signatures.  For  this  purpose,  the  instantaneous  shape 
and  position  of  the  shadows  of  the  surface  signatures  were  projected  on 
a  line-engraved  plexiglas,  placed  at  the  bottom  of  the  test  section,  and 
recorded  on  high-speed  Polaroid  film  through  the  use  of  the  inclined 
mirror.  Each  picture  included  time  accurate  to  0.1  seconds  (generated 
by  an  electronic  timer  and  projected  into  the  inclined  mirror  by  a  small 
magnifying  mirror).  Either  prior  to  or  after  each  experiment,  a  scale 
was  placed  on  the  surface  of  the  water,  and  both  the  scale  and  its 
shadow  were  recorded  on  film.  The  relationship  between  the  actual  scale 
and  its  shadow  was  used  to  transform  the  measurements  obtained  from  the 
shadows  to  those  of  the  actual  disturbances  in  order  to  account  for  the 
parallax,  free  surface  magnification  and  light-source  proximity  effects. 
In  this  manner,  the  actual  separation  of  the  surface  signatures  was 
deduced  from  the  photographs  at  suitable  time  intervals  (determined 
from  the  times  recorded  on  the  films).  All  surface  disturbances  were 
recorded  on  film  until  they  had  completely  dissipated.  The  results  are 
nondi mens ionali zed  and  plotted  in  various  forms  and  compared  with  those 
obtained  theoretically.  Each  experiment  was  repeated  at  least  twice  to 
ascertain  that  the  results  were  reproducible. 


III.  ANALYSIS 


A.  DIMENSIONAL  ANALYSIS 

The  dependent  parameters  of  major  importance  are  the  instantaneous 
positions  of  the  vortex  pair  (x,y)  and  the  surface  disturbances  (s,0). 

They  may  be  expressed  as  functions  of  the  following  parameters: 

x  *  f(t,  U,  dQ,  pQ,  dp/dy,  v,  B,  AR,  a,  g,  rfi)  (1) 

y  =  f(t,  U,  dQ,  p0,  dp/dy,  v,  B,  AR,  a,  g,  rQ)  (2) 

and 

s  =  f(t,  U,  dQ,  pQ,  dp/dy,  v,  B,  AR,  a,  g,  rfi,  )  (3) 

in  which  t  represents  the  time,  U  the  velocity  of  the  model,  dQ  the  initial 
depth  of  the  vortex  pair,  pQ  the  reference  density  of  water,  dp/dy  the 
linear  density  gradient,  v  the  kinematic  viscosity  of  water,  B  the  base 
width  of  the  model,  AR  the  aspect  ratio  of  the  model,  a  the  angle  of 
attack  of  the  model,  g  the  gravitational  acceleration,  a  the  surface 
tension  of  water,  and  rg  an  effective  core  radius,  characterizing  the 
effect  of  the  wing-tip  shape  in  addition  to  other  wing  parameters. 

The  height  and  width  of  the  test  section  and  the  height  of  the  model 
within  it  were  not  included  in  the  foregoing  because  a  detailed  analysis, 
based  on  ideal  vortices,  has  shown  that  the  velocities  induced  by  the 
bottom  were  negligible.  Effects  of  surface  tension  on  the  instantaneous 


position  of  the  vortex  pair  is  deemed  neglible.  It  is  for  this  reason 
that  o  is  not  included  in  Eqs.  (1)  and  (2). 

A  dimensional  analysis  of  Eqs.  (1)  through  (3)  yields 

x/B  =  f (Ut/B,  dQ/B,  NB/U,  U2/gB,  UB/v,  AR,  a,  r0/B)  (4) 

y/B  =  f(Ut/B,  d0/B,  NB/U,  U2/gB,  UB/v,  AR,  a,  re/B)  (5) 

and 

s/B  =  f (Ut/B,  dQ/B,  NB/U,  U2/gB,  UB/v,  AR,  a,  re/B,  P()BU2/a)  (6) 

in  which 

N  *  (-9/P0  dp/dy)1/2  (7) 

and  is  known  as  the  Brunt-Vaisala  frequency. 

Equations  (4),  (5)  and  (6)  may  be  recast  in  terms  of  the  initial 
separation  bQ  of  the  vortex  pair,  and  the  initial  mutual  induction 
velocity  VQ  of  the  vortices  by  noting  that  B/bQ  and  VQ/U  are  uniquely 
determined  by  AR  and  a  for  a  given  wing  shape.  Thus,  one  has 

*/b0  -  f(V0t/b0.  d0/b0>  Nb0/V0,  vjj/gb,,.  V0b0/v,  re/b0)  (8) 

x/b0  -  f(V0t/b0>  d0/b0,  Nb0/V0,  v2/gb0,  V0b0/v,  re/b0)  (9) 

and 


s/bQ  *  f(»0t/b0,  d0/b0.  Nbn/V0,  v2/gbn,  v„b„/v,  r./bn.  Pgb^/o)  (10) 


Experiments  have  shown  that  [Ref.  6]  the  ratio  of  the  initial  vortex- 
core  spacing  to  wing  span,  bQ/B,  is  0.70  for  the  sharp-edged  Delta  wing 
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0.71  for  the  round-edged  Delta  wing,  and  nearly  equal  to  tt/4  for 
the  rectangular  wing  (NACA  0012). 

The  previous  measurements  have  also  shown  that  the  dimensinal  para- 
2 

meters  VQ/gb0  (a  Froude  number  which  ranged  from  0.003  to  0.03)  and 
VQb0/v  (a  Reynolds  number  which  ranged  from  3500  to  10,500)  are  not 
important  within  the  range  of  the  parameters  encountered  in  the  present 
experiments.  With  the  foregoing  arguments  and  noting  that  in  experiments 
in  a  homogeneous  medium,  N=0,  Eqs.  (8),  (9)  and  (10)  may  be  written  as 

»/»0  =  f(V/b0,  d0/b0,  re/b0)  (11) 

y/b0  =  f(V0t/b0.d0/b0,  re/b0)  (12) 

and 

s/b0  =  f<y/b0,  /b0,  re/b0,  P0/b0v2/c)  (13) 

Whereas  the  parameters  VQt/b0  and  dQ/b0  may  be  changed  independently, 
re/bQ  is  taken  as  nature  provides  it.  The  primary  reason  for  this  is 
that  a  century  of  theoretical  and  experimental  aerodynamics  research 
has  been  incapable  of  describing  the  details  of  the  structure  of  the 
tip  vortex  to  be  used  as  initial  conditions  in  a  viscous  solution.  It 
is  surprising,  but  true,  that  until  recently  the  importance  of  the  wing- 
tip  shape  and  its  influence  upon  both  the  initial  tangential  velocity 
profile  and  the  initial  turbulence  in  the  vortex  has  not  been  fully 
appreciated.  Here  the  said  influence  has  been  characterized  in  terms 
of  an  effective  core  radius  with  full  awareness  of  its  shortcomings. 
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In  Eq.  (13),  the  last  term  Pq^qVq/o  represents  the  so-called  Weber 
number.  It  expresses  the  ratio  of  the  inertial  force  to  the  surface 
tension  force  and  hence  the  effect  of  the  surface  tension  on  the  char¬ 
acteristics  of  the  surface  disturbances.  It  is  a  well  known  fact  that 
the  capillary  waves  have  a  minimum  length  of  0.681  inches  (1.73  cm), 
a  minimum  phase  speed  of  9.13  inches/second  (23.2  cm/s),  and  a  minimum 
group  velocity  (at  which  the  energy  is  transmitted)  of  about  7  inches/ 
second  (17.8  cm/s)  [Ref.  17].  The  speed  of  the  surface  disturbances 
observed  in  the  present  investigation  ranged  from  0.5  inch/second 
(1.26  cm/s)  to  4  inches/second  (10  cm/s).  Furthermore,  the  present 
investigation  has  conclusively  shown  that  the  surface  signatures  are 
not  capillary  waves  but  rather  disturbances  generated  by  and  slaved  to 
the  vortices.  Hence  the  effect  of  the  Weber  number  on  the  speed  and 
position  of  the  surface  signatures  is  considered  negligible.  The 
amplitude  of  the  disturbances,  which  may  be  affected  by  the  surface 
tension,  has  not  been  studied  in  the  present  investigation.  The  fore¬ 
going  does  not  rule  out  the  possibility  that  for  larger  vortex  strengths 
the  speed  of  the  surface  signatures  may  coincide  or  exceed  that  of  the 
capi 1 1  ary  waves . 

B.  AN  INVISCID  FLOW  ANALYSIS  OF  THE  MOTION  OF  TRAILING  VORTICES 

It  is  possible  to  analyze  the  interaction  of  a  vortex  pair  with  the 
free  surface  only  when  suitable  assumptions  are  made  regarding  the 
vortices  and  the  free  surface  conditions.  One  can  assume  a  steady  two- 
dimensional  potential  flow  past  a  submerged  line  vortex  and  examine  the 
behavior  of  the  resulting  waves.  Such  an  analysis  has  been  carried  out 
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by  Salvesen  and  von  Kerczek  [18]  for  the  purpose  of  exploring  the  effect 
of  higher  order  approximations  on  nonlinear  problems.  Their  analysis  is 
valid  only  for  waves  with  very  long  lengths  and,  thus,  is  not  of  partic¬ 
ular  interest  for  the  present  study.  Furthermore,  it  is  noted  that  in 
the  Salvesen-von  Kerczek  analysis  the  vortex  is  assumed  to  remain  at 
rest  in  a  steady  stream  of  arbitrary  velocity  U.  It  is  easy  to  show 
that  such  a  situation  is  physically  impossible  and  the  vortex  may  be 
kept  at  rest  only  when  the  free  stream  has  the  proper  velocity  so  as  to 
cancel  the  mutual  induction  velocity  of  the  vortex  and  its  image. 

A  relatively  simple  alternative  is  to  assume  the  free  surface  to  be 
undeformable  at  first  and  calculate  the  resulting  flow  field  due  to  the 
vortex  pair  and  its  image,  and  then  examine  the  behavior  of  the  free 
surface  due  to  the  velocities  induced  on  it.  For  the  case  under  con¬ 
sideration  where  the  disturbances  are  of  very  small  amplitude  and  wave 
length,  the  proposed  method  yields  remarkably  satisfactory  results. 

For  an  inviscid  line-vortex  pair  approaching  the  free  surface 
(assumed  to  be  rigid),  (see  Fig.  2),  the  complex  potential  is  given  by 

w(z)  =  -  j-  [ln(z-zy)  +  ln(z+zv)  -  ln(z-zy)  -  ln(z+zy)]  (14) 

where  r  and  zy  represent,  respectively,  the  strength  and  position  of  the 
vortex. 

The  complex  velocity  anywhere  in  the  field  is  obtained  by  differen¬ 
tiating  Eq.  (14)  with  respect  to  z 

dw  ir  1  1  1  1 

—  =  u-iv  = - (  - + - - - — )  (15) 

dz  2n  z-zy  z+zy  z-zy  z+zy 
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Substituting  z  =  x+iy  and  zy  =  xy+iyv  in  Eq.  (15)  and  separating  the 


real  and  imaginary  parts,  one  has 


<lx  r  y-y  y+y 

dt  2 it  (x+x  +  (y-yur  (x-xjz  +  (y+y„r 


(x+xj2  +  (y+yj2  (x-xj2  +  (y-yj2 


dy  r  x-x  x+x 

v  *  —  *  —  (  - p - *  +  - p - 2 

dt  2 it  (x-x  )z  +  (y-y.,r  (x+xjz  +  (y+y„r 


(x-xy)2  +  (y+yv)2  (x+xy)2  +  (y-yv)2 


since  the  initial  mutual  induction  velocity  is  given  ty 


Eqs.  (16a)  and  (16b)  can  be  recast  as 
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(16b) 


(18a) 


The  velocity  distribution  on  the  free  surface  may  be  found  by  setting 
y/bQ  =  0  and  v/VQ  *  0.  Then,  Eqs.  (18a)  and  (18b)  reduce  to 


where  u$  is  the  velocity  at  the  free  surface. 
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It  is  of  special  interest  to  know  the  position  xm  of  the  maximum 
surface  velocity  u$m  relative  to  the  position  of  the  trailing  vortex. 
The  position  of  the  vortex  is  given  by  [Ref.  19]. 
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<W 


<Vb0)2 


=  4  + 
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(20) 


where  it  is  seen  that  the  relationship  between  xv/bQ  and  yv/bQ  depends 
on  d0/bQ.  Differentiating  Eq.  (19)  with  respect  to  x/bQ,  setting  it 
equal  to  zero,  and,  finally,  combining  it  with  Eq.  (20),  one  has 
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(21) 


where  n  =  xv/b0  and 


k  =  4  +  (b0/d0)2 


The  velocities  at  the  center  of  a  vortex,  located  at  x  *  xy  and 
y  *  yv,  may  be  obtained  from  Eqs.  (18a)  and  (18b)  by  omitting  the  con¬ 
tribution  of  the  vortex  itself  to  the  velocity  at  its  center  [i.e.,  by 
omitting  the  last  term  in  Eq.  (18a)  and  the  first  term  in  Eq.  (18b)]. 
This  procedure  yields 
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(22a) 


It  is  of  particular  interest  to  know  the  normalized  positions  of  the 
vortex  and  of  the  maximum  surface  velocity  as  a  function  of  the  normalized 
time  VQt/b0.  Noting  that 


and  vy 


and  combining  with  Eqs.  (22a)  and  (22b)  one  can  numerically  evaluate  the 
position  of  the  vortex  and  of  the  maximum  velocity  as  a  function  of  the 
normalized  time  VQt/b0  and  the  initial  depth  of  the  vortex  pair  do/bQ. 
This  enables  one  to  determine  first  the  surface  velocity  profile  and 
then  the  position  of  the  maximum  velocity  as  functions  of  VQt/b0,  d0/bQ, 
and  x/bQ.  The  computer  programs  for  the  calculation  of  us/VQ  and  x^bg 
are  given  in  Appendices  A  and  B.  The  program  for  the  simultaneous  cal¬ 
culation  of  xv/bQ,  yv/bQ,  xm/b0»  and  un/Vo  is  Presented  in  Apprendix  C. 
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The  position  of  a  marker  particle  initially  located  at  (0,  dQ)  may  be 
calculated  as  a  function  of  time  for  various  values  of  d0/bQ.  Noting 
that  for  a  fluid  particle  on  the  vertical  axis  x  =  0  and  u  =  0,  Eq.  (18b) 
may  be  written  as 


Writing  v  *  dy/dt  and  noting  that  y  =  dQ  at  t  *  0,  a  numerical  integration 
scheme  can  be  employed  to  calculate  the  position  of  the  marker  particle 
through  the  use  of  Eqs.  (22)  and  (23)  as  a  function  of  dQ/b0  and  VQt/b0. 
The  computer  code  presented  in  Appendix  D  enables  one  to  perform  the  said 
calculations.  The  results  obtained  in  this  manner  are  discussed  later. 

Similarly,  one  can  calculate  the  position  of  the  top  of  the  recircula¬ 
tion  cell  as  a  function  of  VQt/b0  and  dQ/b0  by  noting  that  the  top  of  the 
recirculation  cell  is  a  stagnation  point.  The  code  for  such  a  calculation 
is  given  in  Appendix  E. 


IV.  RESULTS 


A.  PHYSICAL  CHARACTERISTICS  OF  THE  SURFACE  DISTURBANCES 

The  physical  appearance  of  the  surface  signatures  produced  by  a  vortex 
pair  approaching  the  free  surface  can  be  characterized  in  two  distinct 
modes:  "striations"  and  "scars". 

The  striations  appear  abruptly  sometime  after  the  passage  of  the  lift¬ 
ing  surface  as  ridges  perpendicular  to  the  path  of  the  model  (see  (Fig.  3). 
They  grow  rapidly  to  a  total  length  of  about  2.6  bQ  to  3.0  bQ,  (L/2/b0  = 

1.3  to  1.5).  The  pure  striation  phase  of  the  surface  condition  lasts  for 
a  time  interval  of  (VQt/b0)  =  0.6  to  1.0.  Shortly  thereafter,  two  distinct 
signatures  or  'scars'  appear.  They  are  nearly  parallel  to  the  model  track 
and  run  along  the  tips  of  the  striations  (see  Figs.  4  and  5).  While  the 
scars  become  the  dominant  features  of  the  free  surface,  the  striations  . 
remain  in  the  signature  band  and  stretch  more  or  less  randomly  in  the 
transverse  direction. 

The  scars  move  side  ways,  i.e.,  the  width  of  the  signature  band  increases 
with  time,  at  a  more  or  less  constant  speed  if  unimpeded  by  other  phenomena 
such  as  Crow  instability  or  vortex  breakdown  (see  Figs.  6  and  7).  The 
lateral  motion  of  the  scars  is  nearly  in  phase  with  that  of  the  vortices. 

As  the  scars  and  striations  age,  they  attenuate,  spread  over  a  larger  area, 
and  eventually  become  indistinguishable.  The  foregoing  takes  place  at 
relatively  small  depths  (say  for  dQ/bQ  less  than  about  4).  For  larger 
depths  of  submergence,  the  Crow  instability  sets  in  and  amplifies  rapidly 
during  the  later  stages  of  the  rise  of  vortices.  In  the  absence  of  ambient 
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turbulence,  as  in  the  present  case,  the  Crow  instability  is  a  consequence 
of  the  Helmholtz  and  Rayleigh  instabilities  which  lead  to  the  production 
of  turbulence  around  the  periphery  of  the  vortex  core  and  in  the  recircula¬ 
tion  cell.  The  proximity  of  the  free  surface  may  impede  the  occurrence 
of  the  Crow  instability  due  to  the  mutal  interaction  of  the  vortex  pair 
with  its  image,  i.e.,  the  mutual  interaction  of  four  vortices.  Leaving 
aside,  for  the  time  being,  the  exact  nature  of  the  Crow  instability  under 
the  influence  of  the  free  surface,  it  is  evident  from  Figs.  8  through  10 
that  the  cause  of  the  scars  is  in  fact  the  vortices  beneath  the  free 
surface.  In  other  words,  the  scars  and  striations  are  not  small  amplitude 
surface  waves  propagating  independently  away  from  their  generator  but  are 
local  disturbances  whose  generation,  growth  and  propagation  are  controlled 
by  the  vortices.  Figures  8  through  10  are  representative  examples  of  the 
surface  signatures  which  could  have  resulted  only  from  the  vortices  which 
have  undergone  Crow  instability.  The  simultaneous  observations  of  the 
vortex  core  (i.e.,  the  dye)  and  the  surface  signatures  have  shown  that  the 
signatures  are  trapped  by  and  slaved  to  the  vortices  beneath  the  surface. 

Figure  9  shows  a  succession  of  ring-shaped  scars  generated  following 
the  linking  of  the  vortices.  It  is  noted  that  while  the  inside  of  these 
relatively  crude  rings  are  covered  with  randomly  distributed  striations, 
the  region  outside  the  scars  is  disturbance  free.  Figure  10  shows  scars 
similar  to  those  discussed  above  for  a  different  test  run.  It  is  clear 
from  Fig.  9  that  the  linking  does  not  necessarily  occur  simultaneously 
along  the  trail  at  equal  Intervals.  As  was  shown  in  previous  investiga¬ 
tions  [Ref.  6],  the  linking,  like  the  Crow  instability,  is  dependent  on 
the  ambient  turbulence  and  one  cannot  expect  the  linking  to  take  place 
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with  perfect  regularity.  Figures  8  through  10  show  the  nearly  perfect 
symmetry  of  the  ring-shaped  scars  with  respect  to  the  axis  of  the  channel. 
These  figures  also  show  a  number  of  dark  circular  regions  corresponding 
to  local  surface  depressions.  They  result  from  the  simple  fact  that 
when  a  vortex  ring  approaches  the  free  surface  it  breaks  and  opens 
up  like  a  horseshoe  or  a  U-shaped  magnet  under  the  influence  of  its  image. 
Then  the  U-shaped  vortex  immediately  joins  its  mirror  image  at  the  free 
surface,  forming  a  new  ring  (half  real,  half  imaginary),  in  accordance 
with  the  well-known  vortex  theorems.  The  local  surface  depressions 
correspond  to  the  points  where  the  real  and  image  vortex  arms  shake 
hands  at  the  free  surface.  Figure  11  shows  the  scar  and  the  striations 
generated  by  an  inclined  vortex  ring  just  prior  to  the  occurrence  of 
the  breaking-relinking  phenomenon. 

B.  EXPERIMENTAL  DATA 

The  experimental  data  evaluated  in  the  manner  described  previously 
have  been  plotted  in  terms  of  the  governing  parameters.  The  relative 
horizontal  half  spacing  s/bQ  is  shown  as  a  function  of  the  normalized 
time  VQt/b0  in  Figs.  12  through  16  for  the  Round-edged  Delta  wing  (Delta-2R), 
in  Figs.  17  through  20  for  the  sharp-edged  Delta  wing  (Delta-2S,  and  in 
Figs.  21  through  25  for  the  rectangular  wing  (RP-2),  for  the  dQ/bQ  ratios 
noted  in  each  figure.  Each  symbol  in  each  figure  represents  one  inde¬ 
pendent  test  run.  At  least  two  test  runs  were  performed  for  a  given  set 
of  initial  conditions. 

The  characteristics  of  the  data  will  be  discussed  through  the  use  of 
two  representative  plots.  Figure  14  shows  the  data  obtained  with  Delta-2R 
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for  d0/bQ  *  4.  This  figure  reveals  the  following  facts:  (i)  the  surface 
signatures  in  the  form  of  striations  do  not  come  into  existence  until  a 
definite  time  of  VQt/b0,  ( VQt/bo  =  3.3  in  Fig.  14);  (ii)  for  small  values 
of  VQt/b0,  corresponding  to  the  striation  phase,  the  data  show  larger 
scatter;  (iii)  the  pure  striation  phase  of  the  surface  signatures  ends 

(as  determined  from  the  pictures)  after  a  time  interval  of 

AVot/bo  =  °*6 

to  1.0,  (in  Fig  14,  AVQt/b0  =  0.9);  (iv)  there  is  a  particular  time 
interval  where  the  data  represent  the  scar-dominated  surface  signatures 

(from  VQt/bo  *  4.3  to  7  in  Fig.  14);  and,  (v)  the  scatter  in  the  latter 

region  is  not  due  to  the  shift  of  the  data  in  the  vertical  direction 
from  one  run  to  the  next  but  rather  due  to  the  random  errors  as  noted 
from  the  figure. 

Data  for  a  larger  dQ/b0  are  shown  in  Fig.  16.  In  addition  to  all  of 
the  characteristics  noted  earlier,  one  observes  that  the  data  in  Fig.  16 
show  larger  scatter.  This  is  primarily  due  to  the  occurrence  of  the  Crow 
instability  and  the  formation  of  the  vortex  rings.  They  give  rise  to 
sinusoidal  and  ring-shaped  scars  and  render  the  accurate  evaluation  of 
s/bQ  somewhat  more  difficult. 

In  the  foregoing,  all  the  data  obtained  with  the  three  models  at 
various  depths  have  been  presented  (the  tabulated  data  are  given  in 
Appendix  F).  It  has  been  stated  on  the  basis  of  observations  and  mea¬ 
surements  that  the  surface  signatures  are  slaved  to  the  vortices. 

However,  the  fundamental  question  of  what  gives  rise  to  the  disturbances 
and  how  they  are  related  to  the  velocity  field  generated  by  the  vortices 
has  not  yet  been  discussed.  It  has  been  surmised  that  the  position  and 
the  amplitude  of  the  disturbances  are  directly  related  to  the  strain 
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field  arising  from  the  velocity  field  generated  by  the  vortices.  It  is 
with  this  hypothesis  in  mind  that  the  velocity  distributions  obtained 
from  the  inviscid  flow  analysis  will  be  presented  and  subsequently 
related  to  the  experimental  data. 

C.  RESULTS  OF  THE  INVISCID  FLOW  ANALYSIS 

Figures  26  through  31  show  the  surface  velocity  profiles  at  various 
times  for  various  depths,  as  obtained  through  the  use  of  Eq.  (19). 
Evidently,  the  maximum  velocity  becomes  nearly  constant  for  x/bQ  larger 
than  about  2.  This  correponds  to  the  state  where  the  vortex  and  its 
image  are  no  longer  affected  by  the  other  vortex  of  the  pair  and  its 
image.  In  other  words,  the  surface  velocity  for  x/bQ  larger  than  about 
2  is  dictated  primarily  by  one  of  the  trailing  vortices  and  its  image. 

It  is  exactly  for  the  same  reason  that  the  shape  of  the  velocity  profile 
becomes  increasingly  more  symmetrical  about  a  vertical  axis.  Further¬ 
more,  it  is  noted  for  future  reference  that  the  half-width  of  the  velocity 
profile  at  40  percent  of  its  maximum  is  C/bQ  =  0.66. 

Figure  32  is  a  plot  fo  the  relative  position  of  the  maximum  surface 
velocity  as  a  function  of  the  normalized  time  for  various  dQ/b0  values, 
as  obtained  through  the  use  of  Eq.  (21).  For  a  given  d0/bQ,  the  maximum 
surface  velocity  occurs  first  at  a  distance  xm  which  is  larger  than  the 
initial  vortex  separation  (e.g.,  at  xm/bQ  =  1.2  for  dQ/b0  =2).  As  time 
increases  and  as  the  vortex  pair  rises,  x^bjj  decreases.  Subsequently, 
xm/bQ  approaches  a  minimum  and  then  increases,  more  or  less  uniformly. 

The  foregoing  results  are  predicated  on  the  assumption  that  the  flow  is 
inviscid  and  there  is  no  mutual  induction  instability. 


Figures  33  through  38  show  xm/bQ  and  um/VQ  as  functions  of  VQt/b0 
for  various  values  of  d0/bQ.  Also  shown  in  these  figures  (with  dashed 
lines)  are  the  variations  of  xv/bQ  and  yv/bQ  with  the  normalized  time. 
These  figures  reveal  several  interesting  facts  which  will  be  pointed 
out  with  reference  to  Fig.  33.  First,  it  is  noted  that  the  vortices 
are  closer  to  the  vertical  axis  than  the  position  of  the  maximum  velocity, 
i.e.,  xy>xm.  This  rather  surprising  condition  prevails  for  a  time  period 
dependent  on  d0/bQ,  as  evidences  from  Figs.  34  through  38.  Second,  beyond 
a  minimum  v0t/bQ,  the  position  of  the  maximum  velocity  and  the  x-position 
of  the  vortex  are  indistinguishable.  In  other  words,  the  maximum  velocity 
arises  only  from  the  contributions  of  the  vortex  and  its  image  directly 
above  it.  Finally,  the  amplitude  of  the  surface  velocity  increases 
rapdily  with  time  and  reaches  a  nearly  constant  value  when  the  x-positions 
of  the  vortex  and  the  maximum  velocity  become  nearly  identical.  Had  the 
vortices  migrated  from  an  infinite  depth  towards  the  free  surface,  they 
could  not  have  come  a  distance  closer  than  bQ/2  to  the  free  surface  and 
the  surface  velocity  would  have  been  u$  =  2(r/2Tib0/2))  =  4VQ.  In  Fig.  33, 
un/Vo  asymptotically  reaches  a  value  slightly  larger  than  4.0  because 
the  vortices  are  generated  at  a  relatively  small  depth.  In  fact,  Eq.  (20) 
for  xy/b0  *  reduces  to 

yv  2 

(24) 
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which  shows  that  yv/bQ/2  is  smaller  than  1.0  for  vortices  migrating 
upwards  from  a  finite  depth  of  dQ.  For  larger  values  of  d0/bQ  (see  e.g., 
Fig.  38)  um/VQ  asymptotically  approaches  4.0. 

Figure  39  shows  the  position  of  the  marker  particle,  originally  located 
at  (x,y,t)  =  (0,do,0),  relative  to  the  instantaneous  position  of  the 
vortices.  Note  that  at  an  arbitrary  time  t,  the  marker  particle  is 
located  at  (0,y)  and  the  vortices  at  (xy,yv).  The  relative  position  of 
the  marker  particle,  as  given  by  (yv-y)/bQ,  increases  rapidly  at  first 
and  then  reaches  a  maximum  value  equal  to  the  semi -mi nor  axis  of  the  re¬ 
circulation  cell  (see  insert  in  Fig.  39).  For  a  vortex  pair  migrating 
in  an  infinite  fluid,  the  semi -major  and  semi -mi nor  axes  are  given  by 
m/bQ  *  1.045  and  n/bQ  =  0.865,  respectively  [Ref.  19].  For  vortices 
migrating  in  a  finite  depth  of  fluid,  the  smaller  the  initial  depth  of 
the  vortices,  the  smaller  is  n/bQ  (see  Fig.  39).  Accordingly,  m/bQ 
increases  since  the  recirculation  cell  contains  the  same  body  of  fluid 
(the  effects  of  entrainment  and  detrainment  through  recirculation  cell 
boundary  will  not  be  discussed  here).  Once  the  marker  particle  reaches 
the  top  of  the  recirculation  cell,  it  continues  to  move  with  the  recir¬ 
culation  cell,  provided  the  latter  is  not  deformed  due  to  the  free 
surface  effects.  Under  these  circumstances  (yv-y)/bQ  remains  constant 
for  a  certain  period  of  time  (as  seen  from  Fig.  39,  particularly  for 
larger  values  of  dQ/b0).  As  the  recirculation  cell  begins  to  'see'  the 
free  surface,  it  flattens  against  the  free  surface  and  thus  (yv-y)/bQ 
decreases.  Ultimately,  the  vortices  reach  their  maximum  elevation,  as 
dictated  by  Eq.  (21),  and  the  marker  particle  nearly  touches  the  free 
surface,  as  does  the  top  of  the  recirculation  cell. 
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So  far  the  discussion  has  centered  around  the  motion  of  a  two- 
dimensional  vortex  pair.  The  trailing  vortices  are  generated  at  a 
speed  of  U  (model  velocity)  and  rise,  at  least  initially,  with  a  mutual 
induction  velocity  VQ.  Thus  the  vortex  pair,  as  viewed  from  the  side, 
is  inclined  at  an  angle  of  y  =  arctan(VQ/U)  relative  to  the  horizontal. 
Consequently,  the  top  of  the  recirculation  cell  does  not  arrive  at  the 
free  surface  at  all  points  along  the  trailing  vortex  at  the  same  instant. 
In  fact,  if  one  were  to  consider  an  imaginary  point  depicting  the  inter¬ 
section  of  the  top  of  the  recirculation  cell  with  the  free  surface,  one 
would  see  it  move  at  a  speed  of  U  in  the  direction  of  motion  of  the  model 

Figure  40  shows  the  position  of  the  marker  particle  relative  to  the 
initial  position  of  the  vortex  pair  for  various  dQ/bo  values.  In  other 
words,  the  marker  particle  located  at  (0,dQ)  at  t  =  0  moves  to  a  point 
(0,y)  at  time  t.  Thus,  h  =  dQ  -  y  represents  the  migration  of  the  marker 
particle  from  its  original  position.  Since  the  marker  particle  very 
quickly  reaches  the  top  of  the  recirculation  cell.  Fig.  40  shows,  for 
all  intents  and  purposes,  the  position  of  the  top  of  the  cell  relative 
to  the  initial  position  of  the  vortex  pair. 

D.  ANALYSIS  AND  DATA 

It  has  been  hypothesized  on  the  basis  of  the  observations  and  measure 
ments  that  the  creation  and  propagation  of  the  striations  and  scars  are 
a  consequence  of  the  strain  field  arising  from  the  nonuniform  surface 
velocity  distribution.  The  following  is  a  critical  examination  of  this 
hypothesis  through  the  use  of  the  experimental  data  and  the  inviscid  flow 
analysis. 


The  variables  of  primary  concern  in  the  comparison  of  the  theoretical 
and  experimental  results  are  the  time  of  creation,  the  instantaneous 
shape  and  position,  the  rate  of  propagation,  and  the  duration  of  the 
surface  signatures. 

The  data  obtained  with  all  the  models  have  shown  that  the  striations 
always  appear  first  when  the  vortex  pair  arrives  at  a  particular  depth 
from  the  free  surface.  The  measurements  and  the  calculations  based  on 
Eq.  (26)  and  Eqs.  (14)  through  (23)  have  shown  that  the  said  distance  is 
yv/bQ  *  1.0.  The  pure  striation  phase  lasts  a  time  period  of  AVQt/b0  = 
0.6  to  1.0.  During  this  time  interval,  the  amplitude  of  the  surface 
velocity  increases  rapidly  and  nearly  attains  its  maximum  value  (see 
Figs.  33  through  38).  This  rapid  transition  period  culminates  with  the 
appearance  of  scars  as  shown  in  Figs.  4  and  5. 

Figures  41  and  42  are  two  representative  plots  showing  the  position 
of  the  scar  front,  as  deduced  from  the  data,  and  the  position  of  the 
maximum  surface  velocity,  as  given  by  Eq.  (21),  as  a  function  of  the 
normalized  time.  The  conclusions  to  be  arrived  at  on  the  basis  of  these 
two  figures  are  uniformly  valid  for  all  the  data  obtained  during  this 
investigation. 

It  is  noted  from  Figs.  41  and  42  that  (i)  the  variation  of  xsf/bo 
with  VQt/b0  is  essentially  linear,  i.e.,  the  scar  front  has  a  fairly  con¬ 
stant  velocity,  (0.43  VQ  in  Fig.  41  and  0.53  VQ  in  Fig.  42);  (i i )  the 
position  of  the  maximum  velocity  (based  on  the  inviscid  flow  analysis) 
trails  the  scar  front  for  a  normalized  time  smaller  than  about  5.8  (in 
Fig.  41)  and  leads  the  scar  front  for  larger  times;  (iii)  the  scars  come 
Into  existence  at  a  time  when  the  vortex  pair  almost  reaches  its  maximum 
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height  and  begins  to  move  essentially  in  the  lateral  direction  (see  also 
Figs.  33  through  38);  and  finally,  (iv)  the  scars  cease  to  exist  when 
the  vortices  completely  dissipate.  The  normalized  time  of  dissipation 
has  been  previously  established  [Refs.  6,  15  and  16]. 

The  conclusions  reached  on  the  basis  of  Figs.  41,  42,  and  all  others 
like  them,  point  out  one  important  fact;  If  the  position  of  the  scars 
is  slaved  to  the  vortices,  as  shown  by  numerous  test  runs,  and  if  the 
vortices  are  observed  to  dissipate  in  a  finite  time,  then  an  analysis 
based  on  constant  vortex  strength  and  inviscid  flow  assumption  cannot 
adequately  predict  either  the  magnitude  or  the  location  of  the  maximum 
surface  velocity  or  the  positions  of  the  vortices  and  the  scars.  Conse¬ 
quently,  one  must  consider  a  circulation  model  which  accounts  for  the 
dissipation  of  the  vortex  partly  due  to  the  mixing  of  the  counter-sign 
vorticity  in  the  overlapping  regions  of  the  vortex  pair  in  the  recircula¬ 
tion  cell,  partly  due  to  turbulence,  and  partly  due  to  the.  diffusion 
of  vorticity  generated  beneath  the  free  surface. 

Squire  [Ref.  20]  suggested  that  a  turbulent  trailing  vortex  could  be 
described  by  Lamb's  solution  [Ref.  19]  if  the  kinematic  viscosity  v  in 
it  were  replaced  by  an  eddy  diffusivity  vt  as 

-=- =  (l-exp(-n2/4a  ))  (25) 

o 

2  2 

where  n  -  r  /trQ,  r  the  radial  distance  from  the  vortex  core  and  aQ  = 
vt/r0  which  is  a  function  of  the  vortex  Reynolds  number,  rQ/v.  For  the 

4 

experiments  reported  herein,  the  vortex  Reynolds  number  is  about  3x10 
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which  corresponds  to  an  aQ  value  of  2x10  ,  according  to  the  data  compiled 

and  reported  by  Govindaraju  and  Saffman  [Ref.  21].  Thus,  Eq.  (25)  reduces 
to 

r 

_  =  !  .  exp(-A/V  t/b  )  (26) 

o 

2 

where  A  =  20(re/bo)  . 

During  the  early  stages  of  the  motion  of  the  vortex  pair,  the  vortex 
core  relaminarizes  while  the  periphery  of  the  core  becomes  increasingly 
more  turbulent.  Neither  the  distribution  of  vorticity  nor  the  amount  of 
circulation  at  a  given  r  is  measured  or  calculated  to  any  degree  of 
accuracy  in  spite  of  a  century  of  aerodynamics  research.  Notwithstanding 
this  fact,  one  can  obtain  some  idea  about  the  magnitude  of  A  as  follows. 

As  the  vortices  begin  to  rise,  the  distribution  of  vorticity  may  be 
assumed  to  be  Gaussian  with  an  effective  core  radius  smaller  than  bQ/2. 

It  is  noted  that  the  vortices  are  still  a  distance  bQ  apart  at  this  time. 
Thus,  it  is  reasonable  to  assume  rg  *  0.4bQ,  (i.e.,  A  =  3.2),  for  the 
vortices  which  have  not  yet  had  time  to  diffuse  much  as  the  time  they 
come  near  the  free  surface  (a)  distance  of  about  bQ/2,  see  Eq.  (24). 

This  is  particularly  true  for  vortices  generated  in  shallower  depths, 
(small  dQ/b0  values).  The  majority  of  the  dissipation  for  these  vortices 
takes  place  when  they  begin  to  move  sideways  during  which  time  they  give 
rise  to  relatively  strong  scars  and  suffer  from  increased  dissipation  due 
to  free  surface  proximity  and  core  distortion. 

The  vortices  migrating  from  larger  depths  have  sufficient  time  for  the 
spread  of  turbulence  in  their  core,  leading  to  a  larger  effective  core 
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radius.  Noting  that  the  top  of  the  recirculation  cell  is  bounded  by  the 
free  surface,  the  effective  core  radius  is  at  least  equal  to  bQ/2 .  How¬ 
ever,  as  the  recirculation  cell  flattens  against  the  free  surface  and  as 
the  vortex  core  becomes  more  elliptical  in  shape,  the  effective  core 
radius  may  acquire  even  larger  values  (say,  rg  =  0.6bQ  for  which  A  =  7.2). 
The  majority  of  the  dissipation  for  these  vortices  takes  place  before 
they  come  near  the  free  surface.  Thus,  they  move  sideways  at  a  slower 
speed,  dissipate  at  a  slower  rate,  and  generate  relatively  weaker  scars. 

It  is  important  to  note  that  no  simple  expression,  such  as  Eq.  (25), 
based  on  the  replacement  of  the  molecular  diffusivity  by  an  approximate 
turbulence  diffusivity,  is  likely  to  account  for  the  entire  spectrum  of 
the  complex  phenomena  ranging  from  turbulent  diffusion  to  Crow  instability 
to  vortex  linking.  Thus,  Eqs.  (25)  and  (26)  are  not  proposed  here  as 
definitive  solutions,  but  rather  as  means  of  exploring  the  observed 
behavior  of  the  vortices  and  its  consequences.  If  the  hypotheses  made 
prove  to  be  workable  on  the  basis  of  these  simple  turbulence  models, 
then  one  can  explore  in  greater  detail  the  phenomena  controlling  the 
mutual  interaction  of  the  vortices  with  the  homogeneous  and  density 
stratified  media  and  with  the  free  surface. 

With  the  objectives  cited  above,  Eq.  (26)  has  been  incorporated  into 
Eqs.  (14)  through  (23)  and  the  resulting  expressions  were  used  to  plot 
Fig.  43  for  various  values  of  A.  It  is  seen  that  the  variation  of  A 
from  3  to  ®  does  not  materially  affect  x^q  for  normalized  times  smaller 
than  that  corresponding  to  the  minimum  x^bg.  The  reason  for  this  is 
that  during  the  early  stages  of  the  rise  of  the  vortices,  the  dissipation 
of  vorticity  is  not  significant,  particularly  for  vortices  migrating  in 


a  shallow  medium  as  in  Fig.  43.  For  normalized  times  larger  than  that 
corresponding  to  the  minimum  xm/bQ,  the  slope  of  each  curve  depends 
strongly  on  A.  The  reason  for  this  is  that  the  vortices  come  under  the 
influence  of  the  free  surface  and  dissipate  faster  as  noted  earlier. 

The  slope  of  each  curve  represents  the  lateral  velocity  of  the  point 
at  which  the  maximum  velocity  occurs.  Thus,  it  is  clear  that  the  velocity 
of  the  said  point  may  be  rendered  identical  to  that  of  the  scars  by  choos¬ 
ing  judiciously  the  value  of  A. 

Figures  44  through  46  show  the  effect  of  the  variation  of  A  on  xm/b0, 
um/V0,  xv/bQ,  and  yv/bQ.  These  figures  may  be  compared  with  Fig.  35, 
which  corresponds  to  the  case  of  A  *  ®,  ( i . e . ,  r  =  constant).  It  is  noted 
that  um/VQ  is  relatively  smaller  and  begins  to  decrease  after  reaching  a 
maximum  (see  Fig.  46)  instead  of  asymptotically  rising  to  a  constant  value 
(see  Fig.  35). 

One  is  now  in  a  position  to  explore  the  consequences  of  the  hypothesis 
that  the  scars  are  slaved  to  the  vortices  and  that  the  strain  field  arising 
from  the  nonuniform  velocity  distribution  is  the  cause  of  the  scars.  With 
this  objective  in  mind,  the  representative  data  obtained  with  Delta-2R  for 
d0/b0  3  4.0  are  compared  with  the  corresponding  analytical  results  in 
Fig.  43  and  it  is  found  that  the  mean  slope  of  the  experimental  data 
(representing  only  the  scars)  matches  reasonably  well  the  slope  of  one 
of  the  analytical  curves  for  which  A  3  6.  Figure  47  shows  the  comparison 
of  the  said  slopes.  Even  though  the  ability  of  the  remainder  of  the  data 
to  yield  similar  correspondence  remains  yet  to  be  seen,  it  is  at  once 
clear  from  Fig.  47  that  the  scar  and  the  maximum-velocity  point  are  in 
phase  but  are  not  spatially  coincident.  A  simple  calculation  shows  that 


A  =  (xsf~xm)/b0 ,  i.e.,  the  vertical  distance  between  the  dashed  line 
and  the  mean  line  going  through  the  data,  is  about  0.66.  In  other  words, 
the  front  of  the  scar  is  not  at  the  point  of  maximum  velocity,  but  at 
the  point  where  the  velocity  is  a  certain  fraction  A  of  the  maximum.  From 
the  velocity  profiles  similar  to  those  shown  in  Fig.  28,  it  is  determined 
that  the  velocity  is  40  percent  of  the  maximum,  i.e.,  A  =  0.4,  for 
A  =  (xsfxm)/b0  =  °-66- 

The  foregoing  results  have  been  obtained  on  the  basis  of  Figs.  43  and 
47  alone.  It  will  now  be  shown  that  the  values  of  A  and  A  are  nearly  uni¬ 
versal  constants  of  the  surface  velocity  profiles  for  x/bQ  larger  than 
about  1.5.  For  this  purpose,  let  us  consider  Fig.  4a  where  us/VQ  is 
shown  as  a  function  of  x/b„  for  various  values  of  V  t/b„  for  d  /b„  =  4 
and  A  *  6.  This  and  all  other  similar  figures  for  other  relative  depths 
show  that  (i)  the  amplitude  of  the  maximum  velocity  rises  sharply  with 
increasing  normalized  time  and  then  decreases  gradually,  as  evidenced 
by  the  envelope  of  the  maximum  points,  and  (ii)  the  velocity  profiles 
become  almost  symmetrical  about  a  vertical  axis  passing  through  their 
maximum  for  x/bQ  larger  than  about  1.5*. 

A  detailed  analysis  of  those  profiles  which  have  become  symmetrical 
has  shown  that  they  are  also  similar.  In  other  words,  all  symmetrical 
profiles  under  consideration  may  be  reduced  to  a  single  profile  if  the 
velocities  are  normalized  by  um/2  and  the  lateral  distances  5  (half  of 

The  normalized  times  for  x/b  larger  than  about  1.5  also  correspond 
to  the  scar-phase  of  the  surface  disturbances  during  which  the  vortices 
move  laterally,  more  or  less  independent  of  the  influence  of  their  real 
counterpart. 
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< 
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the  width  of  the  velocity  profile  where  u  =  um/2).  From  such  a  nearly 
universal  profile  it  has  been  deduced  that  X  =  5  =  0.66  corresponds  to 
A  =  0.4. 

Special  experiments  have  been  carried  out  in  order  to  determine  the 
relative  positions  of  the  scar  front  and  the  vortex  axis  (which  corresponds 
to  the  position  of  the  maximum  velocity).  Figure  49  shows  the  picture  of 
the  trailing  vortices  and  the  scars.  The  value  of  X  has  been  deduced  from 
such  figures  and  was  found  to  be  equal  to  0.65  +  0.05.  This  result  compares 
exceedingly  well  with  that  cited  above  and  provides  indirect  evidence  of 
the  assumptions  made  regarding  the  relationship  between  the  vortices  and 
the  scars  and  between  the  vortex  strengths  and  Squire's  model.  Obviously, 
additional  evidence  in  the  way  of  analysis  and  measurements  is  needed  to 
strengthen  the  validity  of  all  aspects  of  the  hypothesis. 

In  the  following,  the  entire  scar  data  will  be  compared  with  the  theo¬ 
retical  predictions,  and  then  the  dependence  of  A  on  dQ/bo  and  the  model 
shape  will  be  discussed. 

Figures  50  through  54  show  the  normalized  position  of  the  scar  front, 
xsf/b0,  as  a  function  of  VQt/b0  for  the  round-edged  Delta  model  (Delta-2R) 
for  various  values  of  dQ/b0.  Also  shown  in  these  figures  are  the  theore¬ 
tically  predicted  scar  fronts  and  the  relative  position  of  the  maximum 
velocit”.  It  should  be  emphasized  that  the  data  and  the  analysis  pre¬ 
sented  in  these  figures  pertain  only  to  the  scars.  Similar  data  and 
comparisons  are  presented  in  Figs.  55  through  58  for  the  sharp-edged 
Delta  model  (Delta-2S),  and  in  Figs.  59  through  63  for  the  rectangular 
wing  (RP2). 
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In  comparing  the  relative  positions  of  the  scar  front  with  the  pre¬ 
dictions  of  the  analytical  model,  it  has  been  discovered  that  a  single 
value  of  A  cannot  account  for  the  changes  brought  about  in  the  scar 
velocity  with  the  depth  of  migration  of  the  vortices.  It  has  been  em¬ 
phasized  earlier  that  the  vortices  which  rise  from  relatively  small 
depths  create  stronger  scars,  while  the  vortices  migrating  from  greater 
depths  give  rise  to  weaker  scars.  It  has  also  been  noted  that  the  scars 
generated  by  stronger  vortices  must  move  faster  laterally  because  the 
surface  velocity  at  X  =  0.65  is  necessarily  larger  than  that  generated 
by  weaker  vortices.  Consequently,  one  should  expect  smaller  A  values 
for  larger  scar-front  velocities,  and  vice  versa.  Furthermore,  as  pointed 
out  in  Ref.  [6]  ,  all  of  the  characteristics  of  a  vortex  depend  both  on 
the  circumstances  under  which  it  is  created  and  on  the  conditions  to 
which  it  is  subjected  afterwards.  This  is  particularly  true  for  vortices 
which  propagate  in  a  turbulence  free  environment.  Since  the  roll-up  of 
the  vortices  varies  from  one  model  to  another,  and  since  the  Helmholtz 
and  Rayleigh  instabilities  depend  upon  the  particular  roll-up  process, 
the  initial  turbulence  resulting  from  them  may  vary  significantly.  This 
in  turn,  affects  the  vorticity  distribution  in  the  core  and  the  subsequent 
evolution  of  the  vortex,  including  the  degree  of  its  susceptability  to 
Crow  instability.  In  other  words,  one  does  expect  variations  in  A  from 
model  to  model  for  the  same  d0/bQ. 

Figure  64  shows  the  dependence  of  A  on  d0/bQ  and  the  model  shape.  As 
anticipated  on  the  basis  of  heuristic  reasoning,  A  increases  with  dQ/b0 
and  shows  a  variation  of  +  0.5  about  a  mean  value,  for  a  given  d0/bQ. 
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Some  of  this  variation  could  have  been  attributed  to  the  errors  in  the 
execution  of  the  experiments  and  in  the  evaluation  of  the  experimental 
data,  had  it  not  been  due  to  the  fact  that  the  values  of  A  for  a  given 
model  are  quite  consistent.  It  must  also  be  noted  that  the  values  of 
A  shown  in  Fig.  64  correspond  to  those  calculated  through  the  use  of  an 
effective  core  radius  ranging  from  about  0.4bQ  to  0.6bQ.  Thus,  the  varia¬ 
tion  of  A  from  model  to  model  for  a  given  dQ/b0  may  be  due  to  the  varia¬ 
tions  in  the  core  radius.  Experiments  with  additional  models  will  have 
to  be  carried  out  to  determine  the  dependence  of  A  on  the  model  shape 
or  the  roll-up  process.  Furthermore,  experiments  in  an  environment 
with  prescribed  ambient  turbulence  may  shed  further  light  on  the  variation 
of  A  with  dQ/b0  and  the  model  shape. 


V.  CONCLUSIONS 


The  investigation  reported  herein  warranted  the  following  conclusions: 

1.  The  trailing  vortices  give  rise  to  surface  signatures  in  the  form 
of  striations  and  scars. 

2.  The  striations  come  into  existence  when  the  vortex  pair  is  at  a 
distance  equal  to  one  initial  vortex  separation  from  the  free  surface. 

The  duration  of  the  pure  striation  phase  is  roughly  equal  to  the  time 
that  the  vortex  pair  would  take  to  rise  a  distance  of  about  one-half  to 
one  initial  vortex  separation. 

3.  The  scars  come  into  existence  towards  the  end  of  the  pure  stria¬ 
tion  phase  and  when  the  vortices  are  at  a  distance  of  about  one-half  the 
initial  vortex  separation  from  the  free  surface. 

4.  The  velocity  of  the  scar  front  is  essentially  constant  and 
depends  only  upon  the  initial  position  of  the  vortex  pair  relative  to 
the  free  surface. 

5.  A  theoretical  model  based  partly  on  inviscid  flow  analysis  and 
partly  on  an  approximate  turbulence  model  has  shown  that  the  instantaneous 
position  as  well  as  the  velocity  of  propagation  of  the  scars  may  be  pre¬ 
dicted  with  confidence.  The  single  parameter  expressing  the  relationship 
between  the  eddy  viscosity  and  the  turbulence  prevailing  in  the  vortex 
must  be  obtained  experimentally. 

6.  A  comparison  of  the  characteristics  of  the  surface  signatures 
with  those  of  capillary-gravity  waves  has  shown  that  the  signatures  are 
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Figure  3.  Striations 


Figure  4.  Striations  and  Scars  (Example  1) 
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Figure  5.  Striations  and  Scars  (Example  2) 
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Figure  8.  Crow  Instability  (Example  1) 
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Figure  9.  Crow  Instability  (Example  2) 
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Figure  11.  Vortex  Ring  Approaching  Free  Surface 
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Surface  Disturbance 


[Figure  15.  Surface  Disturbance  Position  for  Delta-2R,  d  /b 


igure  17.  Surface  Disturbance  Position  for  Delta-2S,  d  /b 
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Figure  19.  Surface  Disturbance  Position  for  Delta-2S,  d  /b 
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Figure  21.  Surface  Disturbance  Position  for  RP2,  d  /b 


Figure  25.  Surface  Disturbance  Position  for  RP2,  d  /b 


26.  Surface  Velocity  Profile,  d  /b  =  2,  Inviscid  Model 
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Figure  28.  Surface  Velocity  Profile,  d  ,'b  =  4,  Inviscid  Model 


Figure  31.  Surface  Velocity  Profile,  d  /b  =  7,  Inviscid  Model 


Position  of  Maximum  Surface  Velocity,  Inviscid  Model,  d  /b  =  2  through  7 


e  33.  Position  of  Maximum  Surface  Velocity 
and  Maximum  Velocity,  Inviscid  Mode 
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jjigure  39.  Position  of  Marker  Particle  Inside  Recirculation  Cell 
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Figure  49.  Surface  Scars  and  Vortices 
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Figure  52.  Scar  Front  Position,  Delta-2Rf  d  /b 
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Scar  Front  Position,  Delta-2S, 


jFigure  56.  Scar  Front  Position,  Delta-2S,  d  /b 


Figure  57.  Scar  Front  Position,  Delta-2S,  d  /b 
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Figure  59.  Scar  Front  Position,  Delta-2S,  dQ/b 
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APPENDIX  A:  COMPUTER  CODE  FOR  CALCULATING  u  /V 
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APPENDIX  E:  COMPUTER  CODE  FOR  CALCULATING  POSITION 
OF  TOP  OF  RECIRCULATION  CELL 
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APPENDIX  F:  TABULATED  DATA 
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m  CORRESPONDS  TO  STRIflTION  MODE 
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1.56 

4.01 

1.75 

3.70 

1.51 

3.76 

1.69 

4.96 

2.23 

3.85 

1.65 

4.27 

1.88 

5.27 

2.39 

4.16 

1.81 

4.84 

2.10 

4.58 

2.02 

5.36 

2.33 

5.26 

2.23 

5.83 

2.42 

5.72 

2.60 

6.35 

2.56 

6.55 

2.74 
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I 


1.55 

u- 

3.55 

U- 

s/b0 

Vat/bt 

i  8^o 

VQt/» 

1.01  M 

3.67 

0.78m 

3.4! 

1.12m 

3.73 

1.02m 

3.5! 

1.80 

5.18 

1.85 

4.5 

1.99 

6.32 

2.33 

5.01 

2.22 

6.73 

2.48 

5.4: 

2.33 

5.71 

2.53 

6.21 

6.5! 

0.25m 

0.97  m 

1.46m 

1.78 

1.90 

2.09 

2.28 

2.38 


3.67 

0.65  m 

4.09 

0.75m 

3.43 

0.49m 

4.04 

0.75m 

4.65 

1.50  m 

3.89 

1.02  m 

4.55 

1.35  m 

5.28 

1.70 

4.30 

1.51 

5.28 

1.65 

5.69 

1.90 

5.16 

1.70 

5.69 

1.85 

6.16 

2.15 

5.63 

2.02 

6.06 

2.00 

6.88 

2.30 

6.51 

6.97 

2.48 

2.67 
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DELT fl-2R 


B-  5.02* 


a- 12° 

u-  3.51 
V0t/bQ  «/bQ 

3.74  0.90m 
4.51  1.50m 

4.97  1.75 
5.78  2.20 
5.38  1.90 
6.82  2.60 
7.12  2.70 


dQ/bQm  4.0 

V0/U-  0.0450 

U-  3.59 

V0t/b0  a/bp 

3.72  0.70m 
4.03  1.14m 
5.13  1.74 
5.82  2.06 
6.76  2.40 


b0-  3.700 

U-  3.59 
V0t/b0  »/bp 

3.46  0.83m 

3.93  1.12m 
4.30  1.50m 

4.93  1.56 
5.45  1.90 
5.87  2.16 
6.18  2.21 
7.02  2.40 
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—~m  s  ^-fc  « 


DELTR-2R 

do/bo" 

5.0 

B- 

5.02  * 

ot- 12 

0 

v0/u- 

0.0450 

bo‘ 

3.700 

u- 

3.55 

U- 

3.55 

U* 

3.55 

V/b< 

Z  8/,^o 

V0t/b0  s/b^ 

V/b= 

e/bQ 

4.53 

0.76  m 

5.05 

0.81" 

4.74 

0.89" 

5.09 

1.11" 

5.11 

0.85" 

5.05 

1.16" 

6.65 

1.83 

5.93 

1.49" 

6.25 

1.78 

7.06 

2.05 

6.40 

1.56 

6.36 

1.87 

7.48 

2.22 

6.82 

1.82 

7.28 

2.05 

7.89 

2.40 

7.33 

2.06 

8.51 

2.56 

7.39 

2.27 

8.21 

2.40 

u- 

3.55 

u- 

3.51 

u- 

3.55 

V/b0 

i  "7bQ 

V/bfl 

.  »/b0 

V'b9 

»/bQ 

4.83 

0.92" 

4.78 

0.66" 

4.87 

0.83" 

5.09 

0.96" 

5.29 

1.17" 

5.13 

0.99" 

5.87 

1.27" 

5.85 

1.46" 

5.60 

1.71 

6.38 

1.50" 

6.01 

1.50" 

5.91 

1.74 

6.96 

1.71 

6.78 

1.88 

6.42 

1.85 

7.11 

1.77 

6.82 

1.88 

6.99 

1.96 

7.53 

1.80 

7.55 

2.02 

7.15 

2.20 

7.98 

2.42 

7.56 

2.26 

8.34 

2.48 
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i 


i 


B-  5.02 


B-  5.02 


DELTR-2R  dQ/b0-  6.0 


o 

1 

*— » 
to 

0 

vQ/u- 

0.0450 

i 

o 

_Q 

3.70 

U- 

3.55 

u- 

3.55 

u- 

3.63 

V0t/b0  s/bQ 

VQt/b0  s/bg 

V/b9 

«/bQ 

6.00 

0.60  m 

5.89 

0.77m 

5.83 

0.90  m 

7.70 

1.65 

6.30 

1.07  m 

6.10 

1  .08  m 

8.47 

2.05 

6.71 

1.41m 

7.02 

1.80 

7.28 

1.58 

7.68 

2.02 

7.49 

1.80 

8.35 

2.47 

7.90 

1.84 

8.61 

2.56 

8.27 

1.97 

U-  3.62 

V0t/b0  s/bQ 

5.58  0.67  m 
6.04  0.81m 
6.86  1.93 
7.27  2.16 
8.09  2.52 
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0ELTR-2R  d0/bG-  2-°  B-  5.02  * 

o-10°  V0/U-  0.0400  bQ-  3.700 


u- 

3.67 

U- 

V/bo 

s/b0 

VQt/b, 

1.62 

0.84  m 

1.27 

1.81 

1.05  m 

1.54 

3.05 

1.83 

2.45 

4.09 

2.43 

2.63 

3.81 
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5.50  U-  3.53 

s/bp  V0t/bQ  »/bQ 

0.69m  1.51  0.68n 

1.05m  1.65  0.89m 

1.67  2.79  1.73 

1.80  2.98  1.80 

2.48  3.89  2.29 

4.17  2.35 

STRIRTION  MODE 


DELTR-2R 


d0/b0-  3.0 


B-  5.02 


o- 10 


U-  3.51 


2.82 

3.19 
4.05 
4.28 
4.78 
5.60 

6.19 


s/b 


v0/u-  0.0400 

U-  3.51 


3.01 

3.23 

4.14 

4.37 

5.42 

5.56 


a/b, 


bQ-  3.71 

U-  3.50 
V0t/b0  a/b 

2.59 

2.72 
3.18 
3.54 
4.27 
5.13 

5.72 


U-  3.60 
V0t/b0  a/bc 


U-  4.00 
V0t/b0  a/bc 


2.46 

0.68m 

2.48 

2.79 

1.24  m 

3.13 

3.06 

1.55 

3.88 

3.43 

1.60 

4.44 

3.89 

1.79 

5.00 

4.26 

1.86 

5.79 

4.77 

2.06 

5.18 

2.25 

5.64 

2.41 

0.88m 

2.34 

0.52m 

1.35  m 

2.65 

1.11m 

1.55 

3.17 

1.59 

2.00 

3.68 

1.70 

2.35 

3.89 

2.00 

2.47 

4.31 

2.19 

4.72 

2.41 

5.24 

2.57 

5.76 

2.68 

6.33 

2.94 
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DELTR-2R  d0/b0-  5-°  B-  5.02* 

o-10°  V0/U-  0.0400  bQ-  3.700* 


U-  3.59 


U-  3.52 


U-  3.55 


Vfll^i/ba 

5.82  0.99m 
6.01  1.27m 

6.-47  1.82 
6.94  2.19 
8.15  2.54 
8.71  2.93 


Vnt/bn  «/bQ 

4.89  0.66m 
5.16  0.99m 


5.89 

6.26 


1.77 

1.88 


6.99  2.04 
6.62  1.93 


7.40 

7.72 


2.26 

2.54 


V_t/b„  «/b. 


5.11 

5.11 

5.48 

6.22 

6.49 
6.95 
7.46 
7.93 


0.66  m 
0.83m 
1  .38  m 
1.77 
1.99 
2.13 
2.32 
2.62 


8.70  2.76 
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B-  5.02 


DELTR-2R  d0/b0-  6.0 


a- 10 

o 

vQ/u- 

0.0400 

bo- 

3.70 

U- 

3.52 

U- 

3.55 

u- 

3.55 

V0t/b0  s/ba 

Vflt/b0  ./fa* 

VflVb 

8/bQ 

6.31 

0.77  m 

7.18 

1  .39m 

6.36 

1.28  m 

6.68 

1  .09  m 

7.55 

1.50m 

6.82 

1.60 

7.51 

1.43  m 

7.83 

1.66 

7.28 

1.87 

7.92 

1.55 

8.24 

1.73 

7.74 

1.92 

8.20 

1.67 

9.12 

2.01 

8.20 

2.12 

8.98 

1.82 

9.58 

2.14 

8.66 

2.27 

8.61 

1.76 

9.95 

2.32 

9.03 

2.41 

9.40 

1.87 

10.32  2.46 

9.49 

2.63 

10.41  2.25 
11.01  2.41 


U-  3.55 
VQt/bQ  «/bQ 

6.49  0.73« 

6.86  0.98m 
7.46  1.47  m 
7.92  1.92 
8.06  2.01 
8.43  2.25 
8.89  2.44 
9.26  2.46 
9.63  2.67 
10.22  2.83 
10.64  2.97 
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1.74  0.68m 
1.78  0.86m 

2.17  1.57 
2.65  1.63 

3.18  1.86 
3.48  2.10 
3.83  2.28 
4.13  2.33 
4.87  2.63 


1.17  0.63m 

1.61  1.17m 

1.87  1.47  m 

2.31  1.71 

2.74  1.77 
2.91  1.83 

3.18  2.01 
3.52  2.25 
4.22  2.53 


1.30  0.87  m 

1.74  1.39m 
1.91  1.59 

2.52  1.62 
2.78  1.83 
3.13  2.08 
3.48  2.13 
3.83  2.37 
4.09  2.55 
4.44  2.61 
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DELTA-2S  d0/b0-  3.0  B-  5.02 ' 

a- 10°  V0/U-  0.0330  bQ-  3.700 


u- 

4.06 

u- 

4.06 

v0t/b0 

s/bQ 

V/b0 

«/bQ 

2.61 

0.82  m 

2.35 

0.68m 

3.00 

1.57 

2.52 

0.04m 

3.44 

1.63 

3.57 

1.68 

3.87 

1.77 

3.31 

1.51 

4.05 

1.98 

4.39 

1.98 

4.74 

2.28 

4.52 

2.04 

5.22 

2.45 

5.26 

2.39 

5.83 

2.57 

5.79 

2.45 

6.05 

2.69 

6.26 

2.69 

6.70 

2.81 
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B-  5.02' 


DELTB-2S  d0/b0-  4.0 

o-10°  VQ/U-  0.0330  bQ-  3.700 


U-  4.06 
Vnl/bQ  «/bQ 

3.18  0.78* 
3.35  0.98* 
3.78  1.31* 
4.05  1.37* 
4.57  1.70 
4.83  1.76 
5.52  2.23 
5.87  2.37 
6.13  2.48 
6.48  2.65 


U-  4.06 


3.44 

3.78 

4.00 

4.35 

4.65 

4.96 
5.22 
5.83 
6.09 
6.70 

6.96 
7.31 


•/br 
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DELTA-2S  d0/bQ-  5.0  B-  5.02 ' 

0.-IO0  V0/U-  0.0330  bQ-  3.700' 
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RP2 

<Vb0- 

2.0 

B- 

4.50' 

a- 12 

0 

VQ/U- 

0.0285 

bo* 

3.53 

u- 

3.63 

u- 

3.67 

U- 

3.55 

v/b< 

if£ a 

V/bs 

V0t/bQ 

e/bQ 

1.12 

0.74* 

1.38 

0.71* 

1.45 

0.60* 

1.37 

0.87* 

1.74 

1.19* 

1.78 

1.12* 

1.62 

1.38* 

2.10 

1.64 

2.05 

1 .37* 

1.90 

1.58 

2.45 

1.77 

2.68 

1.60 

2.18 

1.64 

2.81 

1.83 

3.01 

1.99 

2.53 

1.77 

3.09 

2.03 

3.55 

2.09 

2.88 

2.03 

3.52 

2.22 

3.97 

2.36 

2.95 

2.21 

3.80 

2.41 

3.23 

2.39 

4.26 

2.48 

3.48 

2.48 

4.62 

2.60 

4.00 

2.67 

3.97 

2.67 

U-  3.55 


U-  3.55 


0.97 

1.25 
1.75 
2.20 
2.63 
3.04 

3.25 
3.61 
3.91 
4.14 


0.72* 

0.94* 

1.36* 

1.50* 

1.79 

1.98 

2.12 

2.39 

2.54 

2.65 


V0t/b0 

1.26 

1.50 

1.79 

2.03 

2.27 

2.54 

2.81 

3.08 

3.29 

3.56 


s/b, 
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u- 

3.63 

U- 

3.61 

u- 

3.63 

'gt/b5 

V/b5 

a/bQ 

4.71 

0.73m 

4.51 

0.79m 

3.86 

0.80m 

4.85 

0.79m 

4.79 

0.91m 

4.41 

1.08m 

5.13 

1.03  m 

5.07 

0.97m 

4.96 

1.12m 

5.41 

1.44  m 

5.35 

1.03m 

5.11 

1.36m 

5.69 

1.50  m 

5.56 

1.62 

5.39 

1.71 

6.15 

1.79 

5.83 

1.85 

6.19 

2.01 

6.29 

1.79 

6.11 

1.97 

6.70 

2.16 

6.53 

1.91 

6.60 

2.08 

6.82 

1.97 

6.88 

2.26 

7.10 

2.20 

7.16 

2.32 

7.31 

2.26 

7.44 

2.50 

7.66 

2.44 

7.65 

2.61 

7.87 

2.50 

7.93 

2.67 

8.21 

2.79 

U-  3.60 


U-  3.60 


4.28 
4.34 
4.73 
5.04 
5.22 
5.41 
5.77 
5.86 

6.28 
6.40 
6.83 
7.04 


0.60  m 
0.65m 
0.95m 
1.02  m 
1.48m 
1.51 
1.79 
1.86 
1.99 
2.09 
2.21 
2.31 


VflHSs 

4.13 

4.56 

4.71 

4.92 

5.32 

5.52 

5.83 
5.95 
6.38 
6.68 

6.83 
6.89 


0.70m 

0.78m 

1.02m 

1.10m 

1.47m 

1.51 

1.69 

1.81 

1.93 

1.95 

2.26 

2.26 
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RP2 

o-12° 


B-  4.50  - 


d0/b0-  6.0 
V0/U-  0.0285  b0-  3.534  ' 


U-  3.99 
V0t/b0  »/bQ 

5.08  0.77 n 
5.29  1.03  m 
5.99  1.11m 
6.33  1.79 
6.13  1.92 
6.73  2.06 
6.97  2.16 
7.21  2.25 
7.51  2.31 
7.81  2.18 
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